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We  present  an  approach  toward  improving  the  output  signal  to  noise 
power  ratio  of  a  discrete  time  matched  filter  by  employing  a  continuous 
time  filter.  The  approach  exploits  the  bandlimited  nature  of  the  unknown 
ccntinuous  tire  signal,  and  results  in  a  filter  which  is  robust  to  this 
inexact  knowledge.  We  furthermore  show  that  the  robustness  property  does 
not  severely  compromise  performance;  in  fact,  the  resulting  filter  possess¬ 
es  an  output  signal  to  noise  power  ratio  which  upper  bounds  that  of  the 
corresponding  discrete  time  filter. 

I .  INTRODUCTION 


The  matched  filter  has  been  of  practical  interest  ior  some  time.  This 
filter,  which  maximizes  the  output  signal  to  noise  power  ratio,  requires 
knowledge  of  the  signal  input  for  its  design.  In  many  cases  it  is  reason¬ 
able  to  expect  that  the  signal  will  be  known  at  various  discrete  instants, 
thus  admitting  the  design  of  the  discrete  time  filter.  Unfortunately, 
while  it  is  reasonable  to  expect  chat  the  input  signal  will  be  known  at 
discrete  instants,  it  is  another  matter  co  assume  that  it  will  be  known 
exactly  as  a  closed  form  analytical  expression  over,  for  example,  an  inter¬ 
val  of  time,  which  would  be  necessary  for  the  design  of  a  matched  filter  in 
continuous  time.  Design  of  the  continuous  time  filter  is  thus  inhibited  by 
such  inexact  knowledge  of  the  signal  input. 

While  the  signal  may  be  incompletely  known,  it  is  reasonable  to  expect 
that  in  many  cases  it  could  be  modeled  as  bandlimited.  If  we  fuithormore 
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assume  the  signal  is  known  at  a  fixed  number  of  instants,  we  might  hope 
that  a  continuous  time  filter  could  be  designed  which  is  insensitive  to  the 
remaining  inexactness  in  our  knowledge  of  the  signal.  We  also  might  hope 

that  the  resultant  filter's  performance  would  be  good  when  compared  to  that 
of  the  discrete  time  filter. 

In  this  paper  we  show  that  an  approach  realizing  these  aspirations  is 
possible.  We  show  that  the  resulting  filter  is  extremely  robust  to  inexact 
knowledge  of  the  continuous  time  signal,  and  simultaneously  possesses  an 
output  signal  to  noise  power  ratio  which  upper  bounds  that  of  the  discrete 
time  filter.  These  results  show  that  in  many  cases  it  is  desirable  to  con¬ 
duct  the  filtering  operation  in  continuous  time,  even  though  the  signal 
may  be  known  imperfectly. 


II.  DEVELOPMENT 


Suppose  the  continuous  time  signal  s(t)  is  of  finite  energy  and  possess¬ 
es  a  Fourier  transform  which  vanishes  outside  the  interval  [-W,  W],  where 
W  >  0  is  known,  and  suppose  also  that  s(t^)  =  s^  for  i  =  1,  2,  ...,  n.  We 
will  employ  the  filter  whose  output  at  time  T  is  Jv( t)h(T-t)dt,  where  y(*) 
is  the  received  waveform  and  h(*)  is  the  Green's  function  of  the  linear 
filter  (all  integrals,  unless  indicated  otherwise,  are  taken  over  the  en¬ 
tire  real  line).  For  the  purposes  of  this  discussion  we  will  take  T  to  be 
fixed  and  satisfy  T  >  t^  for  all  i  =  1,  2,  ...,  n.  The  output  signal  to 
noise  power  ratio  then  becomes  r  « 


R(s,  h)  = 


[ f s(t)h(T-t)dt]: 
fh  -u)h (T-t )h (T-u)dt  du 


where  R(*)  is  the  autocorrelation  of  the  wide  sense  stationary  noise.  Note 
that  the  finite  energy  condition  on  s(*)  allows  us  to  place  conditions  on 
h(*)  independent  of  s(*)  so  that  R(*,  •)  is  well  defined.  We  will  consider 
the  class  H  of  all  finite  energy  h(’)  such  that  0  <  fj R(t-u)h(T-t)h(T-u) 
dt  du  <  00 .  Let  S  denote  the  class  of  admissable  signals,  i.e.  finite  ener¬ 
gy  signals  possessing  a  Fourier  transform  which  vanishes  outside  [-W,  WJ 
and  equaling  s^  at  t^,  i  =  1,  2,  ...,  n.  We  then  may  interpret  the  filter 
design  problem  as  essentially  a  choice  of  an  element  h(*)  of  H  for  a  given 
element  s(')  of  S. 


In  practice,  of  course,  we  do  not  know  s ( * ) ;  we  only  know  that  it  be¬ 
longs  to  S.  We  must  therefore  accept  a  certain  degree  of  inexactness  of 
knowledge  of  s(*)  in  our  decision  procedure  when  choosing  h(»).  It  would 
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obviously  be  beneficial  to  make  the  choice  in  a  way  which  is  insensitive 
to  such  inexact  knowledge  and  which  leads  to  good  performance.  In  particu¬ 
lar,  we  might  hope  that  we  could  choose  a  filter  with  performance  upper 
bounding  that  of  the  discrete  time  filter  for  any  s(*)  €  S. 

The  problems  imposed  by  the  presence  of  various  types  of  inexact  know¬ 
ledge  have  been  encountered  in  many  other  contexts.  Frequently,  it  has  been 
beneficial  to  employ  approaches  which  are  "robust"  to  the  inexact  knowledge. 
Such  methods  very  often  involve  saddlepoint  criteria,  e.g.  see  [1  -  4]. 

The  general  success  of  this  work  thus  motivates  a  saddlepoint  approach  to 
our  somewhat  different  problem;  we  therefore  seek  s(*)  and  h(*)  correspond¬ 
ing  to  inf  sup  R(s,  h).  Our  first  result  shows  how  the  filter  may  be  de- 
s€S  h eH 

signed  by  this  criterion. 


Theorem  1:  Suppose  det [ (sine (W( *-t  .))*R( • ) ) (t . ) ] .  .  ?  0,  and  suppose  R(*) 

^  J  i  *  J 

has  a  Fourier  transform  that  is  positive  on  [-W,  W].  Then  we  realize  inf 

n  s€S 

sup  R(s,  h)  =  min  max  R(s,  h)  via  h(t)  =  £  c.  sine (W(T-t ,-t) )  and  s(t)  = 
h eH  seS  h eH  i=l  1  1 


n 

[  £  c.  sinc(W(»-t . ) )*R(») ] (t) ,  where  the  constants  c  are  chosen  so  that 
i=l  1  1  1 

sCt^)  =  s^  for  i  =  1,  2,  ...,  n. 


Proof:  For  any  heH,  let  H(t)  =  h(T-t),  and  denote  the  Fourier  transforms 
of  s,  H,  R  by  s,  H,  R  respectively.  Then  an  argument  similar  to  that  of 
[5]  may  be  used  to  show  that  for  a  fixed  s,  sup  R(s,  h)  is  realized  if  and 

~  ~  heH  ~  ~ 

only  if  s(u))  =  R(oj)  H(uj)  for  almost  all  ta.  We  thus  have  H(u))  =  s(io)  R  (w) 

for  almost  all  w  6  [-W,  W].  It  is  routine  to  verify  that  / s(t)H(t)dt  is 

Q  ^ 

independent  of  the  values  H(*)  takes  on  [-W,  W]  ;  moreover  it  is  straight¬ 
forward  to  show  that  JJ R(t-u)H(t)H(u)dt  du  =  J  H(u>)Hf  a)dw,  which 

is  minimized  for  H(*)  as  above  by  choosing  H(w)  =  0  on  f-*.  .  Thus  sup 

heH 

Choos- 


R(s,  h)  is  realized  by  choosing  H(u>) 


-0 


^  (oj)  R(ui) 


-1 


if  ui  6  [-W,  W] 
otherwise 


ing  h(*)  corresponding  to  the  above  H(p)  for  a  given  s(#),  we  thus  must 

attempt  to  choose  s(*)  to  realize  inf  sup  R(s,  h) .  Noting  that  scaling  s 

s€S  heH 

allows  us  to  hold  the  denominator  of  R(s,  h)  constant,  we  therefore  minimize 

W  W 

r- 

n 


/'s(t)H(t)dt-A/7"l!(t)H(u)R(t-u)dt  du  -  £  X , 
J  1=1  1 


s  ( t .  )  =  - 

(21!) 


J  j  j  s<“)s 


-W  -W 


J 
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R_1(w)ej  (U)+W)td(,!  dw  dt - — 


'//// 


;  (w)s (w)R  1  (w)R-1  (w)eJ  (u)C'+^u)R(t-u) 


-W  -W 


1  n  r~  n 

dui  dii)  dt  du  -  -jr~  E  A  |s(w)eJU  i  doj(  where  A  and  {A.}  are  undetermined 
1 i=l  J  1  i=l 

Lagrange  multipliers.  Using  a  variational  approach  applied  to  s(*),  we 

W  W 


minimize  J 


-(e)  =  —  /  /  /  (s+e*6s)(w)  •  (s+e*6s)(w)R  ^ (ta)e'' ^U>+W^tdio  dw  dt 

s  (2tt)  •/.  J 


-W  -W 


du  do) 


we  obtain 


—  ^  j  J  J  J  (s+E*6s) (w)  •  (s+t  »5s)  (u>)R  1(io)R  1  (u))e^Wt+^U^R(t-u) 

-W  -W 

.  n  -  . 

dt  du  -  ~  £  A  f  (s+e*6s)  (u))eju)  ^dio.  Letting  Jg(0)  =  0  for  all  6s, 
1=1  1  W 

w 

ain  [  f  dt  +  -i-r  [f 

(2tt )"J  J 


W 

dw  dt  -  — ^  III  s(u>)R  1(ta)R  1(o))e^a)t+t°U^R(t-u)dw  dt  du - /  /  /  s  (w ) 

(2n)2JJ{VI  (2tt  )^w 


•R  1(u))R_l(io)ej(u,t4<Ju)R(t-u)dw  dt  du  -  E  A.ejwti  =  0  for  almost  all 

^  i=l  1 

U)€ [ -W,  W] .  This  implies  via  a  straightforward  calculation  ^  s(-w)R  ^(a))  - 

—  s(-w)R  1  (o>)  -  4—  £  A.  e^^i  =  0  for  almost  all  w€[-W,  W],  Thus  s(io)  = 
2tt  2n.  .  i 

i=l 

n  ~  _ 

E  c .  R(oj)  e-^  1  •  I,  „  <<*>)  almost  everywhere  for  some  constants 

i=1  1  l"w»  W1 

{c.)n  where  I r  ,,  ,  (•)  is  the  indicator  function  for  the  interval 

i  i=l  [-W,  W]  n 

[ — W,  W].  We  therefore  obtain  s(t)  =  [  E  c,  sine (W( *-t , ) )*R( • ) ] (t )  for 

i=l 

some  constants  {c^}^^. 

In  order  to  establish  that  this  condition  is  sufficient,  we  note  that 

W  W 

a  straightforward  calculation  shows  Jg(r.)  =  J~(0)  +  rJi(O)  +  —j 

W  W  -W-W 

.6s(w)R_1(w)e)(,j4<l,)tdu.  dw  dt  -  y  II  /  I  6s(w)  6s(w}R'1(w)R~1(u>> 

(2v)2JJ  Ju  Jv 


6s (w)  6s(w)R  1(a))R-1(u)> 


I 
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j  (wt+um) 

e  R(t-u)  doj  du)  dt  du  _<  (0)  for  large  positive  A,  in  view  of  the 


hypothesis.  Thus  the  condition  s(t)  =  [  E  c.  sine (W( *-t . ) )*R( • ) ] (t)  for 

i=l  1  1 

some  constants  {c^}^_^  is  f>oth  necessary  and  sufficient.  The  {c^}*?_j  are 

chosen,  of  course,  to  satisfy  the  constraint  s(t..)  =  s^  for  i  =  1,  2,  ..., 

n,  which  is  possible  since  det [ (sine (W(»-t . )*R(*) ) (t . ) ] .  .  ¥  0.  Finally 

^  1  i  *  J 

we  note  that  inf  sup  R(s,  h)  is  realized  by  choosing  h  corresponding  to  the 
s€S  h gH 


above  s,  i.e.  H(w)  is  almost  everywhere  given  by 
H(w) 


%  ci  Vw,  w](w)e"JWti  if  w  E['w»  w] 


0 


otherwise 


Thus  H(t)  =  E  c.  sine (W(t-t . ) ) ,  which  completes  the  proof. 

1=1  1  1 

It  might  be  of  interest  to  determine  just  how  robust  the  filter  of 
Theorem  1  is.  We  do  know  the  filter  is  robust  in  the  minimax  sense,  but  we 
might  in  practice  desire  further  knowledge  of  the  sensitivity  of  the  perfor¬ 
mance  of  the  filter  to  perturbations  in  the  imperfectly  known  signal.  This 
is  addressed  by  the  following  result. 


Theorem  2:  Let  h(*)  be  the  filter  of  Theorem  1.  Then  R(s,  h)  is  independent 
of  s  6  S. 


/■ 


Proof:  Note  that  it  suffices  to  show  /  s (t )h (T-t)dt  is  independent  of  s  €  S. 

W  W 


■Is 


Nov/  for  each  i,  J  s(t)  sinc(W(t-t^)  )dt  - 


s(a))eJ‘Jtda3y*e(-ti' 

-W  -W 


-t  .w+wt)  j 


W 


dw  dt 


J  s  (ta)  e^i10  dw  =  ^  s(ti)  =  ^  si>  Thus,  fs  (t)h(T-t)dt  = 
-W  J 


—  E  c.  s.,  which  establishes  the  desired  result. 
Wi=l  1  1 


QED 


We  therefore  can  see  that  in  view  of  Theorem  2  the  filter  of  Theorem  1  is 
indeed  very  robust  to  inexact  knowledge  of  the  signal  s(*),  as  long  as 
s(*)  €  S.  This  filter  therefore  accomplishes  our  goal  of  desensitizing  the 
performance  of  the  filter  to  inexact  knowledge  of  the  signal.  The  other 
goal  was  to  accomplish  the  desensitization  in  a  manner  which  resulted  in 
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good  performance;  in  particular,  we  would  hope  that  the  performance  would 
be  at  least  as  good  (and  often  better)  than  that  of  the  discrete  time 
filter.  This  is  illustrated  by  the  following  result: 


Theorem  3:  Suppose  R(")  is  differentiable  in  a  neighborhood  of  t ^  for  each 
i.  Suppose  the  output  signal  to  noise  power  ratio  of  the  discrete  time 
matched  filter  exists  and  is  denoted  by  R^,  and  let  h(*)  denote  the 
(continuous  time)  filter  of  Theorem  1.  Then  R(s,  h)  >  Rfl  for  all  s  €  S. 


Proof:  For  each  positive  real  number  At,  let  P.  .  (t)  =  - — 
-  i.  At  At 


I[ti,  :.+At](t) 

Throughout  the  remainder  of  this  proof  assume  At  is  small  enough  so 


that  for  each  i,  [t^,t.+At]  is  contained  in  the  neighborhood  of  t^  for  which 
R(")  is  differentiable.  Note  that  in  view  of  Theorem  2,  it  suffices  to 
show  R(s,  h)  a  particular  s  €  S.  Moreover,  since  h(*)  maximizes 

performance  over  H  for  a  particular  s  6  S,  it  suffices  to  show  that  there 
exists  a  sequence  of  filters  h  e  H  such  that  R(s,  h„)-*-R  ,  where  henceforth 

lx  tv  U  J-J 

s(*)  corresponds  to  h(*)  according  to  Theorem  1.  Let  h.  (t)  =  E  s.P.  ,  (T-t), 

At  l  x,ot 

Since  h^t(*)  €  H  for  all  At  sufficiently  small,  it  therefore  suffices 


tj+At 


to  show  lim  R(s,  h.  )  =  R  .  Now  /s(t)h.  (T-t)dt  =  E  s.  f  s(t)*^—  dt, 

At->0  U  J  At  i=l  1  J 


and  it  thus  follows  that  lim  f s(t)h.  (T-t)dt  =  E  s2.  Moreover,  f fRit-u) 

At-Xl  *  +At  i=l  1  JJ 

n  n  1  /"*  J/*  | 

h.  (T-t)h.  (T-u)dt  du  =  T.  S  s,  s.  /  /  R(t-u)* - rr  dt  du,  and  the  con- 

At  At  i=lj  =  l  1  3  J  J  (At)2 


t .  t . 

i  J 


ditions  of  the  hypothesis  allow  the  use  of  Leibnitz's  rule,  which  implies 


ff  n  n 

lim  II R(t-u)h.  (T-t)h.  (T-u)dt  du  =  E  E  s,  s.  R(t.-t.).  We  therefore 

At*0  JJ  At  1  i = 1  j  =  1  1  J  13 

conclude  lim  R(s,  h,  )-*-R„,  which  establishes  the  desired  result. 

“*>  M  °  QED 

Note  that  the  additional  condition  for  Theorem  3  of  R(*)  differentiable 
in  a  neighborhood  of  t^  for  all  i  is  quite  mild.  The  continuous  time 
filter  in  fact  often  outperforms  the  discrete  time  filter.  For  example, 
if  R(t)  =  sinc(Wt),  where  W  >  W,  n=2,  and  =  1Tm/^  for  SOIm-*  positive 

integer  m,  then  It  is  straightforward  to  show  that  the  continuous  time 
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f liter  outperforms  the  discrete  time  filter  for  all  (s^  ,  s0 )  +  (0,  0) . 

III.  CONCLUSION 


We  have  introduced  a  method  to  improve  the  output  signal  to  noise 
power  ratio  of  a  discrete  time  matched  filter  by  employing  a  continuous 
time  filter.  This  method  exploits  the  bandlimited  nature  of  the  unknown 
continuous  time  signal  and  is  robust  to  this  inexact  knowledge.  We  have 
furthermore  showed  that  the  method  results  in  an  output  signal  to  noise 
power  ratio  which  upper  bounds  that  of  the  discrete  time  filter. 
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